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ABSTRACT 


An investigation into the feasibility of using a 
solid fuel, dual combustion mode ramjet in a high Mach 
number tactical missile was conducted to determine the 
regimes where this propulsion plant might be superior to a 
conventional solid fuel ramjet. The results of the 
study show that at Mach 6, 80,000 feet, with a temperature 
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I. INTRODUCTION 


The purpose of this study was to investigate the 
feasibility of utilizing a solid fuel, dual combustion mode 
ramjet (DMRJ) as the propulsion plant for a high Mach 
number tactical missile. As such, the primary objective 
was to examine the conditions under which the performance 
of a DMRJ might be superior.to that of a conventional solid 
fuel ramjet (SFRJ) at a Mach number of 6.0 using realistic 
propulsion performance parameters. The secondary objective 
was to experimentally examine the feasibility of sustaining 
combustion in a Mach 2.5 flow. 

There are many potential benefits from the development 
of a hypersonic vehicle. Among them are the nation's 
efforts to construct the aerospace plane for hypersonic 
travel through the upper and lower atmosphere for both 
civil aviation as well as space exploration. Benefits to 
the military include powerful tactical/strategic missiles 
and interceptor/fighter aircraft. 

Past tactical missile designs have relied almost 
exclusively on solid rocket motors. Their simplicity of 
design, lightweight construction and ease of handling and 
storage makes them ideally suited for the environment 
aboard ships. Rockets have very high thrust-to-weight and 
can deliver reasonably high specific impulse. They are, 
however, limited to short ranges in order to have suffi¬ 
cient energy during the end-game of the intercept to 


















defeat defensive maneuvers by the target. The Navy’s 


current suite of air-to-air missiles include the AIM-7M 


Sparrow, the AIM-9M Sidewinder and the AIM-54C Phoenix. 


All of these weapons have relatively short range with high 


thrust-to-weights, capable of attaining speeds on the order 


of Mach 1.5 to 2.5 above launch Mach. 


With the advent of very long range detection radars, 
satellite targeting, over-the-horizon data link targeting 


and long range cruise missiles, these weapons are fast 


becoming incapable of adequately defending the fleet in a 
tactical scenario. New weapons like the Tomahawk and 


Harpoon cruise missiles utilize turbojet propulsion that 


has allowed greatly extended ranges. However they were 
designed for slow, non-maneuvering surface targets and are 


themselves subsonic. They are, therefore, practically 
useless against a tactical threat. What is needed is a 


long range, standoff weapon capable of reaching out 30 to 
100 nautical miles with a high terminal velocity giving it 
sufficient energy to defeat the threat. 


In order to achieve long ranges and high velocities. 


specifically terminal velocity, the vehicle must be powered 


throughout the full range of the flight. Rockets use boost 


and sustain propellant grains that burn out prior to 


intercept, generally leaving the missile to ultimately 
coast to the target. A pulsed mAor is one means of 


extending the range of solid rockets. Another viable 


propulsion device for longer ranges is a high speed, 
airbreathing engine. Below Mach 2.5, the gas turbine 
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family of engines is inherently more efficient and will 
outperform all other air breathing power plants. At speeds 
greater than Mach 3.0, however, ram pressure alone is 
sufficient to provide proper compression of the airflow so 
that the compressor actually becomes a liability [Ref. 1: 
p. 38]. In addition, materials limit the maximum tempera¬ 
ture allowed in the turbine section so that as Mach number 
increases, less and less energy can be added to the flow 
before that limit is reached. Thus, thrust falls off 
rapidly on the high speed end for a turbojet. 

Between Mach 3.0 and 5.0, the ramjet becomes more 
efficient than the turbojet. Although ramjets have been 
around for some time, it has been only recently that their 
primary disadvantage has been overcome. Below Mach 0.9, 
ram pressure is insufficient to compress the flow adequate¬ 
ly for efficient combustion [Ref. 2: p. 140]. Thus, a 
ramjet must be brought up to speed before it can be 
employed. The advent of the integral-rocket-ramjet (IRR) 
has made this possible, by employing a solid propellant 
booster to take the ramjet up to its take-over Mach number. 
Figure 1 presents the three power plants discussed thus 
far. 

As the Mach number increases towards Mach 5.0, shock 
compression to subsonic velocities results in large losses 
in stagnation pressure. Her! also, a temperature limit is 
reached for the combustion process due to dissociation of 
reactants in the flow. In addition, since subsonic 
combustion is most efficient when the incoming Mach number 























(a) Solid-fueled rocket 



Figure 1 - Propulsion systems for tactical missiles: 
(a) solid-fueled rocket; (b) integral rocket ramjet; 
(c) turbojet. [Adapted from Ref. 2: p. 139] 




















is on the order of 0.2 to 0.3, the static temperature and 
pressure in the combustor approach material limits. 
Figure 2 delineates the limits of the typical ramjet flight 
envelope. 

In order to achieve hypersonic velocities, scramjet 
technology is the most feasible solution. The decrease in 
shock losses (reducing the combustor flow only to low 
supersonic velocities) more than offsets the increased 
heat addition losses from adding energy to a supersonic 
flow. The lower static temperature associated with 
supersonic flow ensures very little dissociation will 
occur and the need for cooling is reduced. This, in 
conjunction with lower pressures, reduces the structural 
limitation on the combustor itself. Thus, greater energy 
addition is possible, resulting in higher exhaust veloci¬ 
ties and more thrust. 

There are two types of scramjets. The first is a 
conventional ramjet with supersonic flow through the 
combustion chamber. , This is accomplished by designing an 
inlet without a normal shock so that the flow remains 
supersonic. The second type of scramjet is more appro¬ 
priately named a dual-combustion mode ramjet (DMRJ). This 
configuration involves using a conventional ramjet as a gas 
generator that exhausts into a second stream of supersonic 
bypass air. Mixing and combustion of the fuel-rich 
effluent occurs, enabling further heat release. Figure 3 
illustrates both types of scramjets. The primary problem 
with the combustion process is whether or not sufficient 

















Figure 2 - Performance range of ramjets 
[Adapted from Ref. 3] 

Multiple 

Fuel injectors ejectable rocket nozzles 



Oblique Chin inlet 
shocks 


Supersonic combustor 


(a) Integral rocket liquid-fueled scramjet 
Supersonic diffuser Supersonic combustor 


Subsonic 

diffuser 



• Nozzle 


H 


Subsonic combustor 


Fuel injectors 

(b) Liquid-fueled dual-combustor ramjet. 


Figure 3 - Supersonic combustion systems for tactical 
missiles. [Adapted from Ref. 2: p. 142] 
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residence time and mixing exist for sustained combustion, 
given the speed of the flow and the kinetic rates of the 
process. 

Figure 4 illustrates the ideal performance of the 
propulsion devices discussed thus far, and the flight 
regimes in which they have advantages over each other. 

There are several advantages associated with the use 
of solid fuels in both rockets and ramjets. Liquid fueled 
ramjets (LFRJ) require a positive fuel control, a means to 
atomize the fuel droplets and a means of stabilizing the 
flame. These requirements make the the LFRJ inherently 
more complex, as well as heavier, than the SFRJ. 

The SFRJ is "self-throttling 1 *. The amount of mass 
flow taken into the engine and its properties determines 
the fuel flow rate in the combustor. The general expres¬ 
sion for the regression rate for an SFRJ is: 

(eqn. 1) r = k G x pY t z 

where G represents the mass flux entering the combustor, p 
is the chamber pressure, t is the air inlet (static) 
temperature and x, y and z are specific to the fuel 
and generally lie between 0.1 and 0.8. The flame stabili¬ 
zation is accomplished through the geometry of the combus¬ 
tor rather than insertion of vee-gutters or other bluff 
body flameholders. Figure 5 illustrates the flow field in 
the flame stabilization region of an SFRJ combustor. 
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Figure 4 - Theoretical performance envelopes. 
[Adapted from Ref. 2: p. 144] 























Although the SFRJ is "self-throttling”, the throttling 
is not optimum. The positive fuel controller in the LFRJ 
will enable it to achieve higher performance throughout the 
flight envelope than the SFRJ. This situation is changing 
as metallized fuels enable a higher energy density in the 
SFRJ fuel grain (resulting in improved performance) and 
variable air flow geometries are being developed. An SFRJ 
is able to achieve very high performance with a simple 
design. 











II. NATURE OF THE PROBLEM 

A. ANALYTICAL SOLUTION: SOLID FUEL RAMJET 

For the purposes of this study, analytical models had 
to be assumed for both the SFRJ and the DMRJ. For the 
SFRJ the following conditions were employed: 

(a) Design point: M 0 « 6.0; h = 80,000 ft; 1 lbm/sec 
flow rate 

(b) Four-shock inlet with ramp angles of 12, 28 and 43 
deg aes (from horizontal), or an ideal inlet with 
no shock losses 

(c) Inlet operating 10% supercritical 

(d) Subsonic pressure recovery = 0.95 in the inlet 

(e) Combustor inlet Mach =0.2 

(f) Constant area combustor 

(g) One-dimensional flow in the combustor and nozzle 

(h) HTPB grain as fuel (p = 0.0332 lbm/in 3 ) 

(i) Cq = 1.5 in the SFRJ combustor 

(j) Equilibrium, adiabatic combustion 

(k) Nozzle throat choked 

(l) Isentropic exhaust nozzle 

(m) Frozen flow expansion through the exhaust nozzle 
to ambient pressure 

Figure 6 depicts the station numbers in effect for the 
analysis. 
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The continuity equation (conservation of mass), in the 


form below, was used to size the cowl and the combustor. 


(eqn. 2 ) 


Y g e I 

paV—= 

Y + 1 

RT ' 

Y + 1 / 


A/A * 

M 


Oblique and normal shock relationships were used to 
determine the pressure loss up to station 2. With Mach 
number and flow rate specified and stagnation temperature 
known, equation 2 yields the area, A 2 . 

The only known parameters after the combustor were the 
throat Mach number and the exit pressure. To proceed 
across the combustor, conservation of energy and momentum 
are applied. For the energy equation, the Naval Weapons 
Center (NWC) Propellant Evaluation Program, PEPCODE, was 
used for equilibrium, adiabatic combustion to construct 
tables of combustion gas properties and theoretical 
stagnation temperatures at station 4 as a function of the 
stagnation properties entering the combustor and the 
fuel-to-air ratio (f) [Ref. 3]: 

(eqn. 3) PEPCODE - fn(f,T,P) «> T th , C p , R, 7 

Three tables are included in Appendix A. In order to 
determine the properties at station 4, knowing the proper¬ 
ties from PEPCODE, the following equations had to be 
simultaneously satisfied to determine M 4 : 


20 











(eqn. 4) 


T - T 
4 th 2 


1 - 0.35<t> 092 $ s 1 

(eqn. 5) x\ = 

* T 0.3 + 0.353> a92 ® s 1 


where 


(eqn. 5a) 

(D = —f— 

f 

stoich 

(eqn. 6) 


(eqn. 7) 

u = V Y gR t 

(eqn. 8) 

Pi A l -Pe A e- F D =m ' U '- m i U i 

(eqn. 8a) 

F CoPU 

D ~ 2 RTg c 

(eqn. 9 ) 

Y 

P = p(l+ 

RAMJET, an 

interactive computer code 


solve this system of equations, given the output of PEPCODE 
and an initial guess for M 4 as the input data. Each 
iteration was checked with the continuity equation until a 
specified tolerance was met. In this case, that tolerance 
was 0.0001 in mass flow rate. The exhaust velocity, thrust 
and specific fuel consumption (SFC) were then determined 


from: 












(eqn. 10) 


u. = V2C /c «/T Y 


(eqn. 11) 


m « U '~ m o u o 


(eqn. 12) 


SFC = — 
F 


A sample output is included in Appendix B, along with 
the program. 

The grain length is determined by specifying s> for each 
case. Equation 1 determined the regression rate with, 

k = 7.1789xl0“ 5 
X = 0.53 
y = 0.33 
Z = 0.71 

for HTPB. With fuel flow rate specified, the grain 
length, Lg, is determined from, 


(eqn. 13) 


m, = p ji D L r 
f r r p g 


B. ANALYTICAL SOLUTION: DUAL COMBUSTION MODE RAMJET 

As with the SFRJ, a model had to be created to deter¬ 
mine all the properties through the DMRJ. In their 1958 
paper on supersonic combustion, Weber and MacKay [Ref. 5] 
studied a one-dimensional model of a pure scramjet with a 
constant area combustor. Their model used liquid hydrogen 
as the fuel, for which a fixed combustion efficiency of 95% 
was assumed. They accounted for neither momentum losses 
















due to fuel injection in a high speed flow and shock losses 
in the combustor, nor were wall friction losses applied in 
the combustor. Dissociation was addressed briefly. For 
their study, the scramjet outperformed the conventional 
ramjet over the entire range of combustion temperatures 
examined (taken at an altitude of 120,000 ft. Mach 6) as 
long as there was no recombination in the exhaust nozzle 
expansion process (i.e., frozen flow) [Ref. 5: p. 106,110]. 
This stresses the impact on performance that dissociation 
has due to very high static temperature in the ramjet 
combustor. For an isentropic inlet, their scramjet model 
was superior to the conventional ramjet above Mach 7. With 
a wedge-type inlet, that advantage shifted to Mach 5. 
[Ref. 5: p. 102] 

In this study, an attempt was made to take advantage 
of known realistic performance parameters for momentum loss 
in the SFRJ combustor and known subsonic combustion 
efficiency behavior of the solid fuel, HTPB. A one-dimen¬ 
sional, constant area combustor was also employed. Equili¬ 
brium adiabatic combustion through PEPCODE accounted for 
some of the effects of dissociation, though not completely, 
and yielded variable properties. Rather than assume a 
value of combustion efficiency, rj AT was used as a parameter 
in the performance study. 

Figure 7 depicts a schematic representation of the DMRJ 
used in this study with applicable station nurhers. 
Upstream of the "B" path of the mixer was an SFRJ gas 
generator identical to Figure 6. For use with the DMRJ, 

































































however, station 4 has been renumbered as station 3. 
Two separate inlets were assumed; one for the SFRJ gas 
generator and one for bypass air. All the assumptions from 
the previous section for the SFRJ were applicable to the 
gas generator except in this case, the total design mass 
flow rate was 1 lbm/sec. For the bypass portion, the 
following conditions were employed: 

(a) Three-oblique-shock inlet with the identical 
pressure recovery as the SFRJ inlet (prior to the 
normal shock), and an ideal inlet with no losses 

(b) One-dimensional flow 

(c) Mixer inlet Mach = 2.5 

(d) Constant area supersonic mixer 

(e) No mixing losses 

(f) Static pressure at the mixer inlet determines 
exhaust pressure for the gas generator exhaust 
(i.e. SFRJ nozzle sized to match pressures at 
station 3) 

(g) Equilibrium adiabatic combustion 

(h) Frozen flow expansion to atmospheric conditions 

(i) Isentropic nozzle 

The method for solving the SFRJ portion of the DMRJ 
was identical to the previous section, except that the 
exhaust pressure was now equal to P3 a» vice atmospheric. 
RAMJET was again used to determine all the properties at 
station 3B. As before, oblique shock relationships were 
used to determine the properties at station 3A. Equation 2 
was then used to determine the areas at stations 3A and 3B. 











The mixing process presented the most difficulty in 
the analysis. No information was available in the open 
literature to make a reasonable assumption for the mixing 
loss. The approach taken was to assume a value of thermal 
efficiency (r? AT ) • representing losses in the heat release 
process only, to vary the equivalence and bypass ratios 
and to not account for the mixing loss at all. From this, 
a "baseline” performance could be identified to which 
mixing losses could be applied as they become known. 

PEPCODE was again employed to determine the properties 
at the exit of the mixing chamber, given the stagnation 
temperature of the Mach 6.0 flow, an estimate of the 
chamber stagnation pressure (to be verified) and the 
fuel-to-air ratio (f) . In this analysis, a mass averaged 
temperature was determined for station 3, assuming complete 
mixing and no combustion, defined as: 


(eqn. 14) = T 3 * 6 '** 1 '** + m 3B^iB T 3B 

(m M + m 38 )C p3 


Proceeding as before, with eta specified arbitrarily, all 
the properties known at station 3, and station 4 gas 
properties from PEPCODE, all the flow properties at station 
4 were determined by solving equations 2,4,6,7 and 9 
simultaneously. 

MIXER, another interactive code, was written to solve 
this simultaneous set. Input variables included all the 
properties at stations SA^and 3B, the PEPCODE output for 
station 4 and an estimate of the Mach number at station 4. 





























Rather than using the continuity equation as the check, 
this program determines the velocity of the flow at 
station 4 by two means: (1) by equations 4,6 and 7; 
(2) by conservation of momentum, equation 8. Repeated 
estimates of the Mach number are made until the difference 
between the two velocities were on the order of 0.1% . 
This program and a sample output are included as Appen¬ 
dix c. 

Once the properties at station 4 were determined, the 
exhaust velocity, thrust and specific fuel consumption 
were again determined using equations 10 - 12. This 
calculation was performed for each bypass ratio considered 
and each value of thermal efficiency. 













III. ANALYTICAL RESULTS 


A. SOLID FUEL RAMJET 

The performance of the SFRJ is shown in Figure 8. The 
cases considered were: 1) design Mach number of 6.0, 
4-shock inlet; 2) design Mach number of 6.0, isentropic 
inlet; and 3) design Mach number of 4.0, 4-shock inlet. 
All cases were done for an altitude of 80,000 feet. Tables 
of properties for these engines are included in Appendix D. 

The shape of the curves is highly dependent on the 
adiabatic combustion efficiency, » 7 AT . Figure 9 represents 
a plot of equation 5, which is the efficiency behavior of 
HTPB. The minimum at * = 1.0 explains the pronounced dip 
in the performance curve. This dip is present with or 
without losses, regardless of design Mach number. The 
lower design Mach number, however, is more sensitive to 
changes in v • For example, fc_* an overall equivalence 
ratio of 1.3 (r? AT « 0.75), a 1% change in efficiency at 
Mach 4 will cause a 37% greater temperature rise than for 
Mach 6: 


AT 2243 | 1640 | 

* ?AT ~ AT th " 3011 | 4 ” 2200 | 6 

This temperature rise is directly related to the exhaust 

velocity, therefore the thrust. Thus, trends in the 

performance curve are accentuated at the lower Mach number. 
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SFRJ PERFORMANCE 



Figure 8 - Solid Fueled Ramjet performance for design 
altitude of 80,000 ft. 






















The performance curves indicate that specific thrust 
rises with 4> despite the behavior of which is shown 

in Figure 10. The reason for this is two-fold. The 
primary effect is that of decreasing molecular weight. 
This is reflected in higher values of the gas constant, R, 
for higher values of $ (refer to tables in Appendix A) . 
Equation 10 shows that the dependence of exhaust velocity, 
and therefore, thrust, is J R. In addition, increasing 
$ implies an increase in m e . Thus, both of the terms in 
the jet thrust (m e u e ) increase while the ram drag term 
(m 0 u 0 ) remains constant. 

The secondary effect is that of the increasing combus¬ 
tion efficiency of HTPB in the fuel-rich regime (refer to 
Figure 9) . The combination of these two effects appear to 
offset the sharp decrease in theoretical adiabatic combus¬ 
tion temperature, enabling the thrust to increase. 

Inlet losses also have a very large effect on the 
performance of the SFRJ. Figure 8 shows that for the same 
design Mach number, inlet losses reduce the specific thrust 
of the SFRJ by approximately 2 0% throughout the range of 
equivalence ratios. This is due solely to the stagnation 
pressure drop across the inlet system. 

The Mach 4 curve lies above the Mach 6 curve, as 
expected for a ramjet. The ram drag is the primary 
reason for this, as illustrated below for 4> = l.jf. 
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Figure 10 - Equilibrium adiabatic combustion tempera¬ 
ture for HTPB for isentropic inlet (640 psi) and for 
4-shock inlet (120 psi). 
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Although there is a large increase in exhaust velocity 
in going to Mach 6, there is an even greater increase in 
the freestream velocity, which leads to higher thrust at 


Mach 4. 

A second contributor to the higher performance at the 
lower Mach number is the decrease in shock losses due to 
the weaker shocks at Mach 4. 

Dissociation effects are accounted for by PEPCODE. As 
expected for the SFRJ, an increase in design Mach number 
diminishes the effect of adding fuel to the flow. For 
example, for Mg ■ 2.8, the difference in theoretical 

stagnation temperature rise across the combustor going 
from 4> = 0.5 to $> = 1.0 (i.e., doubling the fuel flow) is 
approximately 1300°R. For M 0 =■ 6.0, however, the same 
increase in fuel flow leads to a difference in theoretical 
stagnation temperature of approximately 580°R. This is the 
effect of dissociation on energy addition potential. 

The unique combustion efficiency behavior for the SFRJ 
results in a predicted F/mg versus $ behavior that is 
somewhat different than the behavior that is usually 
reported for ramjets with fixed combustion efficiency. 

Specific thrust will maximize at very rich equivalence 

ratios. 

B. DUAL-MODE RAMJET 

The results from the analytical model support the 

previous theoretical works of both Billig [Ref. 2] and 
Weber & MacKay [Ref. 5]. The performance of the DMRJ 
















exceeds the performance of the conventional SFRJ, within 
the constraints of the conditions applied, by approximately 
20% in the range of anticipated equivalence ratios (0.8 - 


1.1). Each parameter was examined for its impact on the 
calculations and results. However, because no attempt was 
made to account for mixing/shock/wall friction losses in 
the supersonic combustor, these results identify a minimum 
required performance, or a ••baseline”, to which these 
losses must be applied. 

The parameters examined included the combustion 
efficiency (r? AT ) , the bypass ratio 1 , the design Mach number 
and the overall equivalence ratio. The effects of using 
an isentropic inlet were looked at briefly, as well as the 
impact of dissociation in the supersonic combustor. The 
following bypass ratios were used in the analysis: 50/50 
(i.e., 50% primary air, 50% bypass air), 60/40, 75/25 and 
80/20. Tables of properties through the engines are 
included in Appendix E and will be referred to below. 

It must be emphasized that each point represents a 
different design condition and, therefore, a different 
engine. Each engine was re-sized as a parameter was 
changed and then flown at its design point. (The design 
point was Mach 6, 80,000 feet unless otherwise specified.) 

1. Performance Comparison: r? AT 

As with the SFRJ, the performance of the DMRJ is 
highly sensitive to the combustion efficiency specified. 

1 Bypass ratio is defined as the ratio of mass flow 
through tne primary gas generator to the mass flow through 
the bypass duct. 

































The gas generator has the same efficiency characteristics 
as the SFRJ (refer to Figure 9) . The efficiency in the 
supersonic combustor was varied between 0.25 and 0.9. 

Figure 11 is a plot of specific thrust versus the 
overall equivalence ratio for the conventional SFRJ, the 
DMRJ with an efficiency of 90% for all bypass ratios and 
the DMRJ for a combustion efficiency of 25% for two 
bracketing bypass ratios. 

To understand the effect that efficiencies have on 
the DMRJ, it was necessary to compare an isentropic SFRJ to 
an isentropic DMRJ with »? AT = 1.0 for all combustion pro¬ 
cesses. In this case, the engines were isentropic in that 
there were no inlet losses present. For the table below. 
Mg = 6.0, h = 8 0,000 ft. and 4 * TO t = 1,0 same total 

fuel flow). 

SFRJ _ 80/20 _ 50/50 

*SFRJ i.OO 1.25 2.00 

T 3 — 5600 4848 

T 4 5565 5565 5565 

P 4 593 445 291 

M 4 0.33 2.10 1.48 

|P 4 -Pg| 45 193 347 

T 3ave — 5230 4166 

AT — 335 1399 

t 4 5490 3544 4355 

F/m 0 103.6 101.4 97.4 

What this information shows is that at Mach 6, 

80,000 ft., the SFRJ has superior performance to the DMRJ. 

* 

The reason behind this is the stagnation pressure drop due 
just to heat addition. For the SFRJ with subsonic heat 
addition, the drop is very small. It, therefore has the 


mmmm. 

















m 

& 


3 


THRUST VS OVERALL PHI 


LEGEND 

SFRJ 

50/50 [0:9)' 
60/40 (0.9) 
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Figure 11 - Performance comparison for SFRJ, DMRJ with 
rj AT = 0.9 for four bypass ratios and DMRJ with » 7 AT = 
0.25 for two bracketing bypass ratios. (M 0 = 6.0, 
80,000 feet) 
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highest stagnation pressure. For the 50/50 DMRJ, the 
cooling effect of the bypass air is very significant 
(T 3a ve = 4166°R). With considerable unburned fuel and a 
low inlet temperature, much of the heat addition occurs 
in the supersonic flow. Thus, it has the largest pressure 
drop. The 80/20 DMRJ lies between these two extremes. 
This result supports the finding of Weber & MacKay that for 
the same heat addition (both types of engines have the same 
T 4 ) the engine with the highest P 4 will have the highest 
exhaust velocity, which means highest thrust [Ref. 5: 
p. 104]. 

r; AT will affect the stagnation temperature at the 
exit of the combustor, supersonic or subsonic. When 
f 7 AT = 0.9 in the supersonic flow, the DMRJ exhibits a 
signifigant performance advantage at lower values of 
overall equivalence ratio. However, on the fuel-rich end, 
the DMRJ levels off while the SFRJ continues to increase, 
eventually causing a cross-over in performance. This 
behavior is explained in the table below for a DMRJ with a 
bypass ratio of 80/20: 



J SFR 


J 

SFR 

•J 

4* SFRJ 1 * 1 

73 OT? 

73 |277 

73 

rrs 

73 

4* TOT 0.8 

8 0.9 

0 | 1.7 

6 

1.8 

0 

T 3Bth 539 

1 529 

2 | 458 

1 

503 

1 

T 3B 469 

3 462 

6 | 461 

2 

485 

0 

T 4 519 

5 462 

6 I 500 

8 

485 

0 































performance advantage. The third and fourth columns, 
compare the SFRJ and DMRJ at the fuel-rich end ($tot = 1-8) 
where the performance crosses over. 

For the first case, the SFRJ efficiencies are 
identical so that the temperature at the exit of the 
combustors are nearly the same. The only difference was 
that $ = 1.1 lies closer to the peak of the temperature 
profile, shown in Figure 10, than does $ = 0.9, leading to 
a higher theoretical stagnation temperature, hence a higher 
actual stagnation temperature. T 4 is then significantly 
higher due to very efficient burning of the excess fuel in 
the second combustion process. This is the driver behind 
the performance advantage for the DMRJ at the lean end of 
overall equivalence ratio. 

For the second case, the performance is very nearly 
equal, with the SFRJ being slightly superior. In this 
case, the DMRJ primary combustor is very fuel rich at a $ 
of 2.2. This point is well past the peak in theoretical 
stagnation temperature for HTPB. The DMRJ actually 

exhibits a drop in stagnation temperature at the primary 
burner exit from an increase in the fuel flow while the 
SFRJ exit temperature increases. The efficiency behavior 
of HTPB causes the actual exit temperature to increase 
slightly, but it is still less than the exit temperature 
for the SFRJ. The crossover in performance for the SFRJ 
comes from a combination of higher stagnation pressure and 
smaller molecular weight (larger R) that just offsets the 
difference in T 4 . 



















Figure 12 is a plot of DMRJ specific thrust versus 
overall $ for two bypass ratios with rj AT as a parameter. 
This plot shows the increment by which the performance 
varies with efficiency. It appears that for any realistic 
advantage from the dual-mode to exist, combustion efficien¬ 
cies on the order of 70% should be attained (refer to 
Figures 11 & 12) . In the general region of interest 
(i.e., $ < 1.2), 50% efficiency leads to near equal 

performance to the SFRJ accounting for heat addition losses 
but no mixing or shock losses. 

2. Performance Comparison:Bvpass Ratio 


Both Figure 11 and 12 demonstrate the dramatic 
effect of decreasing bypass ratio as efficiency decreases. 
For very high supersonic combustion efficiency, the 
performance lines for each bypass ratio very nearly 
overlay. However, as that efficiency is decreased, the 
performance lines begin to diverge. For the same value of 
overall $, there is a large difference in performance 
between a bypass ratio of 80/20 and one of 50/50. In the 
table below, a closer examination is made for M 0 = 6.0, 
80,000 ft. with inlet losses and with r? AT = 0.5: 


S’TOT 
S’SFRJ 
7AT I RJ 
m-e 


50/50 

0.90 

1.80 

0.90 

0.066 


80/20 

0.88 

1.10 

0.68 

0.065 





























For very nearly the same fuel flow, there is a 
difference in performance. Even though the 50/50 DMRJ has 
a higher stagnation temperature exhausting into the mixer 
(T 3 ) , due to the higher r? AT , T 3ave is significantly lower 
due to the larger amount of cooler bypass air (Tq = 3198). 
This causes T 4 for the 50/50 DMRJ to be lower than the 
80/20 DMRJ, leading to a higher specific thrust for the 
larger bypass ratio for nearly the same amount of fuel 
addition. 

Another effect of bypass ratio is that to obtain 
high values of overall $, larger bypass ratios must be 
used. For example, to obtain an overall equivalence ratio 
of 1.5 would require $SFRJ “ 1.875 for the 80/20 DMRJ. The 
50/50 DMRJ however, would require a primary combustor 
equivalence ratio of 3.0, for which the efficiency data for 
HTPB is not applicable. Thus, for any specified overall 
equivalence ratio, the 50/50 DMRJ will be much more 
fuel-rich than the 80/20 DMRJ. 

One final observation is that although the combus¬ 
tion efficiency in the primary combustor is higher in terms 
of temperature rise for the fuel-rich engines (lower 
bypass ratios), there will be more unburned fuel present in 
the supersonic combustor. For a high r? AT in the supersonic 
combustor, it doesn't appear to matter, and the performance 
of the two cases should be similar. Figure 12 verifies 
this behavior. For the lower r? AT in the supersonic combus¬ 
tor, that excess fuel is burned less efficiently so that a 


























performance difference should exist, 
also depicted in Figure 12. 


This behavior is 


3. Performance Comparison: M 0 

Figure 13 shows how decreasing the design Mach 

number to Mach 4 (80,000 ft.) affects the performance of 

the SFRJ and the DMRJ. The bypass ratio of the DMRJ 

is 80/20 and the combustion efficiency is 90%. Inlet 

losses are included. As in the case of the SFRJ, going to 
a lower Mach number improves the performance of the DMRJ 
for the same reasons: lower shock losses and reduced ram 
drag. Figure 14 includes an 80/20 DMRJ with r? AT = 0.6 to 
show the effect of changing the efficiency and the Mach 
number. 

4. Performance Comparison: Inlet losses 

The effect of using an isentropic inlet on the DMRJ 
with » 7 AT * 0.9 is depicted in Figure 15. Just as the 

performance for the SFRJ increases, the DMRJ realizes a 
gain in specific thrust of approximately 25% at Mach 6, 
80,000 ft. The reason for the improvement, as before, is 
the absence of stagnation pressure losses due to inlet 
shocks. The SFRJ is included in Figure 15 to show that the 
isentropic DMRJ gives the best performance. 

5. Performance Comparison: Dissociation 

One of the reasons for going to supersonic combus¬ 
tion over subsonic is dissociation of the reactant mole¬ 
cules due to the extremely high static temperatures that 
are encountered when a high speed flow is shocked down to 
low subsonic Mach numbers (0.15 to 0.3). Dissociation 




















MACH NUMBER EFFECT 



Figure 13 - Effect of design Mach number on performance 
for SFRJ and 80/20 DMRJ with r,^ = 0.90. 

(4-shock inlet, 80,000 ft.) 





















































prevents complete combustion of the reactants by tying up 


chemical energy in intermediate species. Rather than a 
complete combustion process, where all of the fuel and 
oxidizer react to form CO 2 and H 2 O, species such as OH, 
CO and H are present as well. Their presence reduces the 
potential temperature rise obtained by burning a specified 
amount of fuel and if recombination of the species is 
prevented in the exhaust nozzle, reducing the exhaust 
velocity. This can occur when kinetic rates and residence 
times are such that equilibrium cannot be attained. The 
limiting case is termed "frozen flow". When kinetic rates 
are fast enough compared to residence time in the combust¬ 
or, "shifting equilibrium" exists, where recombination of 
the dissociated species allows more thermal energy to be 
released. 

In using PEPCODE to determine the properties after 
equilibium adiabatic combustion, two approaches can be 
adopted. In this study, stagnation temperature was used 
as an input variable. This means that the gas properties 
are based on stagnation conditions, i.e., M 4 = 0. For the 
SFRJ, static and stagnation temperatures are very nearly 
the same so that it is not a problem. For the DMRJ, 
however, this is not the case. M 4 is supersonic, so that a 
large difference exists between static and stagnation 
properties. A more accurate method would be to use static 
temperature as the input variable in such a way that 
PEPCODE computes properties, based on shifting equili¬ 
brium, for M 4 * 0. This approach, while more accurate, 
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requires significantly more computation. For each data 
point, PEPCODE must be run twice, as well as RAMJET, to 
solve for the properties at the exit of the combustor. 

This effect was examined and the results are 
presented in Figure 16 for an 80/20 DMRJ, with inlet 
losses, Mach 6, 80,000 ft. and r? AT = 0.9. For $ < 1, using 
properties based on static conditions has no significant 
effect on the performance of the DMRJ. For $ > 1.25, i.e., 
fuel rich, properties based on static conditions would 
result in approximately 5% greater specific thrust. Since 
the anticipated range of equivalence ratios is from 0.8 to 
about 1.1, this effect was neglected and the properties 
were computed using stagnation conditions. 




























DISSOCIATION EFFECT 




































































V. CONCLUSIONS AND RECOMMENDATIONS 


From an analytical standpoint, this study showed that a 
solid fuel, dual combustion mode ramjet exhibits areas of 
potential superior performance over a conventional solid 
fuel ramjet for use in a Mach 6 tactical missile. Based on 
the specified flight conditions and performance parameters, 
not accounting for any losses due to mixing in the super¬ 
sonic combustor, the DMRJ outperforms the SFRJ at near 
stoichiometric equivalence ratios for » 7 AT > 0.5. Addition¬ 
al results are summarized below: 

(1) It is estimated that with mixing losses, combustion 
efficiency in the supersonic combustor must be at 
least 70% or greater to achieve superior 
performance with the DMRJ. 

(2) For supersonic combustion efficiencies less than 
approximately 60%, performance of the DMRJ is very 
sensitive to the amount of airflow in the bypass 
stream. Greater bypass flow led to higher perfor¬ 
mance for the same overall $. 

(3) The effect of decreasing inlet losses and lower 
Mach numbers were as expected; both leading to 
better performance. 

(4) Dissociation effects proved to be relatively minor 

( = 1%) in the range of $ near stoichiometric 

(i.e., * « 0.8 - 1.1). 

In order to refine the range of required combustion 
efficiencies, follow-on research should focus on the 
actual mixing processes, the combustion kinetics and 

supersonic combustion processes. Experimental testing with 
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the completed DMRJ apparatus will provide insight into 
some of these areas, as well as some validation of the 
analytical model. Further studies should be conducted 
using solid fuels with different combustion efficiency 
characteristics than HTPB in order to more closely examine 
fuel effects on --he system performance. 
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APPENDIX A - PEPCODE Data 
TABLE 1 

PEPCODE DATA FOR EQUILIBRIUM ADIABATIC COMBUSTION OF HTPB 


0.5 

.0367 

1.266 

53.256 

.326 

4746 

1.297 

0.6 

.0440 

1.263 

53.426 

.330 

4938 

1.292 

0.7 

.0514 

1.260 

53.689 

.334 

5088 

1.289 

0.8 

.0587 

1.259 

54.013 

.338 

5203 

1.287 

0.85 

.0624 

1.258 

54.199 

.339 

5251 

1.286 

.9 

.0661 

1.258 

54.415 

.341 

5292 

1.286 

.95 

.0697 

1.258 

54.631 

.343 

5325 

1.285 

1.0 

.0734 

1.258 

54.848 

.344 

5353 

1.285 

1.05 

.0771 

1.258 

55.095 

.346 

5375 

1.285 

1.1 

.0807 

1.258 

55.342 

.347 

5391 

1.285 

1.15 

. 0844 

1.258 

55.605 

.348 

5402 

1.286 

1.2 

.0881 

1.258 

55.883 

.350 

5406 

1.286 

1.25 

.0918 

1.259 

56.176 

.351 

5405 

1.287 

1.3 

.0954 

1.259 

56.470 

.352 

5398 

1.288 

1.4 

.1028 

1.261 

57.119 

.355 

5365 

1.290 

1.5 

.1101 

1.263 

57.814 

.357 

5308 

. 1.293 

1.6 

.1174 

1.265 

58.540 

.359 

5230 

1.296 

1.7 

.1248 

1.268 

59.313 

.361 

5136 

1.300 

1.8 

.1321 

1.271 

60.101 

.363 

5031 

1.304 

1.9 

. 1395 

1.274 

60.873 

.364 

4923 

1.308 

2.0 

. 1468 

1.277 

61.676 

.366 

4809 

1.312 

2.1 

.1541 

1.280 

62.480 

.367 

4694 

1.316 

2.2 

.1615 

1.283 

63.268 

.369 

4581 

1.32 

2.3 

.1688 

1.286 

64.056 

.370 

4467 

1.324 


Remarks: 


1) All quantities in standard British units 

2) P 2 = 120 psi; T = 3200°R 

3) M 0 = 6.0; h = 80,000 ft 

4) Exhaust data for frozen flow expansion 

















TABLE 2 

PEPCODE DATA FOR EQUILIBRIUM ADIABATIC COMBUSTION OF HTPB 


3 


0.5 

.0367 

1.265 

53.102 

.326 

4812 

1.282 

0.6 

.0440 

1.261 

53.194 

.331 

5041 

1.276 

0.7 

.0514 

1.257 

53.364 

.335 

5229 

1.272 

0.8 

. 0587 

1.255 

53.612 

.339 

5378 

1.269 

0.85 

. 0624 

1.255 

53.766 

.341 

5437 

1.268 

0.9 

.0661 

1.254 

53.951 

.342 

5489 

1.268 

0.95 

.0697 

1.254 

54.137 

.344 

5531 

1.267 

1.0 

.0734 

1.254 

54.353 

.345 

5565 

1.267 

1.1 

.0807 

1.254 

54.848 

.348 

5606 

1.267 

1.2 

.0881 

1.255 

55.404 

.351 

5614 

1.268 

1.3 

.0954 

1.256 

56.037 

.353 

5586 

1.270 

1.35 

.0991 

1.257 

56.393 

.355 

5559 

1.271 

1.4 

.1028 

1.258 

56.748 

.356 

5525 

1.273 

1.45 

.1064 

1.259 

57.119 

.357 

5484 

1.274 

1.5 

.1101 

1.261 

57.505 

.358 

5437 

1.276 

1.6 

.1174 

1.263 

58.293 

.359 

5332 

1.280 

1.7 

.1248 

1.266 

59.096 

.361 . 

5216 

1.284 

1.8 

.1321 

1.269 

59.931 

.363 

5094 

1.288 

1.9 

.1395 

1.273 

60.749 

.365 

4971 

1.292 

2.0 

.1468 

1.276 

61.568 

.366 

4848 

1.296 

2.1 

.1541 

1.279 

62.403 

.368 

4725 

1.300 

2.2 

.1615 

1.282 

63.191 

.369 

4606 

1.304 

2.3 

.1688 

1.286 

63.994 

.370 

4487 

1.308 

Remarks: 







1) All 

quantities in standard 

British 

units 


2) P 2 

= 620 psi 

; T = 3200°R 




3) M 0 

= 6.0; h 

= 80,000 

ft; isentropic 

inlet 


4) Exhaust data 

for frozen flow expansion 
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TABLE 3 


PEPCODE DATA FOR EQUILIBRIUM ADIABATIC COMBUSTION OF HTPB 


s> 

f 

y 

R 

_Cp_ 


7a 

C np 

0.6 

.0440 

1.267 

52.876 

.323 

__ 

3946 

1.295 

—pe 

.298 

0.7 

.0514 

1.261 

52.913 

.329 

4218 

1.288 

.304 

0.8 

.0587 

1.257 

53.035 

.334 

4436 

1.282 

.310 

0.9 

. 0661 

1.254 

53.295 

.338 

4595 

1.279 

.314 

1.0 

.0734 

1.253 

53.673 

.341 

4694 

1.277 

.318 

1.1 

. 0807 

1.254 

54.183 

.344 

4735 

1.278 

.320 

1.2 

.0881 

1.255 

54.826 

.347 

4718 

1.280 

.322 

1.3 

. 0954 

1.258 

55.570 

.348 

4648 

1.284 

.323 

1.4 

.1028 

1.261 

56.389 

.350 

4546 

1.289 

.324 

1.5 

.1101 

1.266 

57.357 

.351 

4415 

1.294 

.324 

1.6 

. 1174 

1.269 

58.098 

.352 

4302 

1.299 

.324 

1.7 

.1248 

1.273 

58.956 

.353 

4178 

1.304 

.325 

1.8 

.1321 

1.277 

59.806 

.354 

4052 

1.309 

.325 

1.9 

.1395 

1.281 

60.640 

.355 

3932 

1.314 

.326 

2.0 

.1468 

1.285 

61.498 

.356 

3810 

1.319 

.327 

2.1 

.1541 

1.289 

62.333 

.357 

3690 

1.324 

.327 

2.2 

.1615 

1.294 

63.144 

.358 

3576 

1.329 

.328 


Remarks: 

1) All quantities in standard British units 

2) P 2 = 50 psi; T = 1640°R 

3) M 0 = 4.0; h = 80,000 ft 

4) Exhaust data for frozen flow expansion 






















APPENDIX B - RAMJET 
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PROGRAM RAMJET 


PROGRAM USED TO DETERMINE THE EXIT CONDITIONS OF A RAMJET BURNER ** 
CHAMBER. A LEGEND FOR THE VARIABLES INVOLVED IS OFFERED BELOW. ** 


STATION 0: FREESTREAM CONDITIONS 

STATION 2: INPUT CONDITIONS FOR THE COMBUSTION CHAMBER 
STATION 4: EXIT CONDITIONS OF COMBUSTION CHAMBER 
STATION E: EXHAUST CONDITIONS(POST NOZZLE) 


M: MACH NUMBER 

T: STATIC TEMPERATURE (DEG R) 

TT:TOTAL TEMP (DEG R) 

P: STATIC PRESSURE (LBF/IN2) 

PT:TOTAL PRESSURE (LBF/IN2) 

GAMMA: RATIO OF SPECIFIC HEATS 
R: GAS CONSTANT (BTU/LBM-DEG R) 

CP:SPECIFIC HEAT (BTU/LBM-DEG R) 

F: FUEL TO AIR RATIO 

PHI: EQUIVALENCE RATIO 

MOOT: MASS FLOW RATE (LBM/SEC) 

U: FLUID VELOCITY (FT/SEC) 

GC: PROPORTIONALITY CONSTANT (FT-LBF/LBM-SEC2) 
J: CONVERSION FACTOR - 778 BTU/FT-LBF 
DELTA: CONVERGENCE LIMIT 
D: PORT DIAMETER (IN) 

STEP: MACH NO. INCREMENT 
THRUST: UNINSTALLED THRUST (LBF) 

STHRST: SPECIFIC THRUST (LBF-SEC/LBM) 

SFC: SPECIFIC FUEL CONSUMPTION (LBM/LBF-SEC) 
ETA: THERMAL EFFICIENCY (TEMPERATURE RISE) 








rMW 
























DOUBLE PRECISION PO,P2,P4,PTO,PT2,PT4,TO,T2,T4,TTO,TT2,TT4,M2,M4 
DOUBLE PRECISION CP2,CP4,GAMMA2,GAMMA4,R2,R4,F,H,MDOT0,GC,CD,U2 
DOUBLE PRECISION U4,TEST,DELTA,MD0T4,STEP,DELTT,MO,UE.UO,THRUST 
DOUBLE PRECISION STHRST,D2,D4,PI,GAMMAE,ME,TE,PE,PTE,TTE,ETA 
DOUBLE PRECISION TT4TH,J,CPE,PHI,SFC 
CHARACTER*! ANS, ANSWER 

100 PRINT *,'ALL UNITS ARE STANDARD ENGLISH UNITS' 

C.INITIAL DATA 

DATA MO,PO,TO/6,.404,390/ 

DATA CD, GC, DELTA, STEP/1.5,3-2.174, .0001, .00001/ 

DATA M2,GAMMAO.RO,CP0.D2/.2000,1.40,53.3,.24,1.674/ 

DATA R2,CP2,GAMMA2/53.3,.28,1.33/ 

DATA J,PT2.MDOTO/778,146.72,1.OODO/ 

PE - PO 

PI—ASIN(1.D0)*2.DO 
A2- (PI/4.DO)*D2*D2 
PRINT *,'A2«',A2 


-CONSTANT AREA COMBUSTOR 


-INITIAL FLOW COMPUTATIONS 

PTO-PO*(1.+M0*M0*(GAMMAO-1.)*.5)**(GAMMAO/(GAMMAO-1.)) 
TTO-TO*(1.+MO*MO*(GAMMAO-1.)*.5) 


T2-TT2/(1.+M2*M2*(GAMMA2-1.)*.5) 

P2-PT2/((1.+M2*M2*(GAMMA2-1.)*.5)**(GAMMA2/(GAMMA2-1.))) 
PRINT 5,'P2-',P2 

5 FORMAT(IX,A3,F9.4) 

PRINT 6,'T2-',T2 

6 FORMAT(IX,A3.F9,4) 

C.INITIAL GUESS ON EXIT MACH 

50 PRINT *, 'INPUT ESTIMATE OF M4' 

READ *, M4 

C-INPUT PEPCODE DATA 






ISVtvSyis 


m 






>/■.■ -.-vre »j'"'njf'-w*r 
















PRINT *,'INPUT TT4TH(DEG R)' 

READ *, TT4TH 

PRINT *,'INPUT GAS PROPERTIES (F,R4,CP4,GAMMA4,GAMMAE)' 

READ *, F,R4,CP4,GAMMA4,GAMMAE 
PHI - F/0.0734D0 
DETERMINE ETA FOR HTPB 
IF (PHI .LE. 1) THEN 

ETA - 1.DO - .35D0*PHI**. 92D0 
ELSE 

ETA - .3D0 + .35DO*PHI**.92D0 
END IF 

SOLVE FOR M4 

TT4 - ETA*(TT4TH - TT2) + TT2 
IF (M4 .GE. .999) THEN 

PRINT *, 'FLOW IS CHOKED FOR THIS VALUE OF F' 

GOTO 50 
ENDIF 

CPE - R4*GAMMAE/((GAMMAE-1.DO)*J) 

T4-TT4/<1.+M4*M4*(GAMMA4-1.)*.5) 

U4-M4*SQRT(GC*GAMMA4*R4*T4) 

U2-M2*SQRT(GC*GAMMA2*R2*T2) 

-MOMENTUM EQUATION FOR P4 

P4-P2+(MDOTO*(U2-(1.+F)*U4))/(GC*A2)-CD*P2*U2*U2/(2.*R2*T2*GC) 
PT4-P4*(1.+M4*M4*(GAMMA4-1. ) *. 5 ) ** (GAMMA4/(GAMMA4-1.)) 

MDOT4—MDOTO*(1.+F) 

TEST - MD0T4 - P4*A4*U4/(R4*T4) 

-CHECK GUESS ON M4 
IF (ABS(TEST) .GE. DELTA) THEN 
M4-M4+STEP 
GOTO 10 
END IF 

■SOLVE FOR EXIT CONDITIONS 

ME-SQRT((2.DO/(GAMMAE - 1.D0)*((PT4/PE)**((GAMMAE-1.D0)/GAMMAE) 

► * IDO))) 

UE-SQRT(2.*GC*J*CPE*TT4*(1.-(PE/PT4)**((GAMMAE-1.)/GAMMAE))) 











TE-TT4/(1.D0+ME*ME*(GAMMAE-1.DO)*.5DO) 

DO-MO*SQRT<GC*GAMMA0*R2*T0) 

THRUST-(MDOT4*UE-MDOTO*UO)/GC 

STHRST-THRUST/MDOTO 

SFC - F/STHRST 

PRINT 15,'MO','M2','M4','ME' 

FORMAT(/IX,4A9) 

PRINT 20,MO,M2,M4,ME 
FORMAT(IX,4F9.3/) 

PRINT 25,'TO','T2','T4','TE' 

FORMAT(IX,4A9) 

PRINT 30,T0,T2,T4,TE 
FORMAT(IX,4F9.1/) 

PRINT 35,'TTO','TT2','TT4','TTE' 

FORMAT(IX.4A9) 

PRINT 40,TTO,TT2,TT4,TT4 
FORMAT(IX,4F9.1/) 

PRINT 45,' PO', ' P2', ' P4' , ' PE' 

FORMAT(IX,4A9) 

PRINT 60,P0,P2,P4,PE 
FORMAT(IX,4F9.3/) 

PRINT 65,'PTO','PT2','PT4','PTE' 

FORMAT(IX,4A9) 

PRINT 70,PTO,PT2,PT4,PT4 
F0RMAT(1X,4F9.3/) 

PRINT 75,'UO','U2','U4','UE' 

FORMAT(IX,4A9) 

PRINT 80,U0,U2,U4,UE 
FORMAT(IX,4F9.3/) 

PRINT 85,'F','MDOTO','THRUST','SPTHRUST','SFC' 

FORMAT(IX,5A9) 

PRINT 90,F,MDOTO,THRUST,STHRST,SFC 
FORMAT(IX,4F9.4,F10,6/) 

WRITE (*,*) 'PERFORM ANOTHER RUN THIS FLIGHT CONDITION (Y/N)? 1 
READ (*, ' (A) ') ANS 











IF (ANS .EQ. 'N') GOTO 110 
GOTO 50 

110 WRITE <*,*) 'START NEW PROBLEM (Y/N)?' 
READ (*,'(A)') ANSWER 
IF (ANSWER .EQ. 'N') STOP 
GOTO 100 
END 


SAMPLE OUTPUT: 


run rami fortran 

VS FORTRAN COMPILER ENTERED. 14:06:42 

★★RAMJET** END 3F COMPILATION 1 ****** 
VS FORTRAN COMPILER EXITED. 14:06:45 


EXECUTION BEGINS... 

ALL UNITS ARE STANDARD ENGLISH UNITS 

A2- 2.20090103 

P2» 142.8812 

T2-3177.0289 

INPUT ESTIMATE OF M4 

7 

.25 

INPUT TT4TH(DEG R) 

7 

5353 

INPUT GAS PROPERTIES (F,R4,CP4,GAMMA4,GAMMAE) 
7 

.0734,54.848,.344,1-258,1285 


MO 

6.000 

M2 

0.200 

M4 

0.293 

ME 

4.309 

TO 

390.0 

T2 

3177.0 

T4 

45*8.4 

TE 

1261.3 

TTO 

3198.0 

TT2 

3198.0 

TT4 

4598.7 

TTE 

4598.7 

PO 

0.404 

P2 

142.881 

P4 

130.670 

PE 

0.404 

PTO 

637.867 

PT2 

146.720 

PT4 

137.879 

PTE 

137.879 

UO 

5805.821 

U2 

538.372 

U4 

931.042 

UE 

7287.592 

F 

0.0734 

MDOTO 

1.0000 

THRUST 

62.6804 

SPTHRUST 

62.6804 


SFC 

0.001171 

PERFORM ANOTHER RUN THIS FLIGHT CONDITION (Y/N)7 
START NEW PROBLEM (Y/N)7 
R; T-0.97/1.54 14:07:18 










APPENDIX C - MIXER 


PROGRAM MIXER 

CHARACTER ANS 

REAL EIA,M4,T4TH,TT4TH,TT4,T4,FT4,P4,A4,A3A,A3B,GAMMA4,R4,CP4 
REAL TT3A,Tr3B,TT3BAR,U3A,U3B,U4,MD0T3A,MD0T3B,MDOT4,GC 
REAL FNT, FNIT, FNU, FNA, FNGAM, FNP, P3A,P3B,V4, DELPA, DIFF 
-INPUT ALL THE PERTINENT DATA 
PRINT *, 'INPUT R4,GAMMA4,T4TH' 

READ *,R4 / GAMMA4,T4TH 

PRINT *, 'INPUT MDOT3A / MDOT3B / TT3A,TT3B' 

READ *, MD0T3A,MD0T3B,TT3A / TT3B 
PRINT *,'INPUT P3A,P3B,A3A,A3B' 

READ *,P3A,P3B,A3A,A3B 
PRINT *,'INPUT U3A,U3B' 

READ *,U3A,U3B 

PRINT *, 'INPUT ETA,M4' 

READ *, ETA,M4 

GC = 32.174 

TT4TH ■ PNTT(T4TH,GAMMA4 / M4) 

MDOT4 = MDCT3A+MDCT3B 

TT3BAR =(MDOT3A*TT3A + MD0T3B*TT3B) /MDOT4 
TT4 = (TT4TH - TT3BAR) *ETA + TT3BAR 
T4 = FNT(TT4 / GAMMA4,M4) 

U4 = FNU(GAMMA4,R4,T4,M4) 

A4 = A3A + A3B 

PT4 = (MDOT4*SQRT(TT4) *FNA(GAMMA4,M4) )/(A4*FNGAM(GAMMA4,R4)) 
P4 = FNP(M4,FT4,GAMMA4) 

DELPA = P3A*A3A + P3B*A3B - P4*A4 

V4 = (. 9*GC*DELPA + MDOT3A*U3A + MD0T3B*U3B) /MD0T4 

DIFF = (U4 - V4) 

PRINT 15,TT3BAR 
















.<n 


n 



15 

FORMAT(IX,'TT3BAR 

,F9.3) 



PRINT 20.TT4TH 


m 

20 

FORMAT(IX,'TT4TH 

,F9.3) 

m 


PRINT 25.TT4 






Vj- 
■> \ 

25 

FORMAT(IX,'TT4 

,F9.3) 



PRINT 30,T4 


& 

30 

FORMAT(IX,'T4 

—',F9.3) 

mm 


PRINT 35 ,U4 


Wit 

35 

FORMAT(IX,'U4 

= ' ,F9.3) 

Ik 


PRINT 40.PT4 



40 

FORMAT(IX,'PT4 

,F9.3) 


PRINT 45,P4 


P 

r / 

45 

FORMAT(IX,'P4 

,F9.3) 

50 

PRINT 50,V4 

FORMAT (IX, 'V4 

,F9.3) 

i 

55 

PRINT 55.DIFF 

FORMAT(IX,'U4-V4 

,F9.3) 


IF (ABS(DIFF) ,LT. .01) STOP 
PRINT *,'CONTINUE? (Y/N)' 
READ (*, ' (A) ' ) ANS 
IF (ANS .EQ. 'N') STOP 
GOTO 10 
END 


rk1cirirfrk1cickir/rkyc1ck1cJrirtctcJcJrjc*-ii-irA-A-&-k-jc-)c 


*** THIS FUNCTION DETERMINES STAGNATION TEMPERATURE (DEG R) GIVEN *** 
*** STATIC TEMPERATURE, MACH AND GAMMA. *** 



FUNCTION FNTT(T,GAM,M) 

£ 

REAL T.GAM.M 


FNTT- T*(l.+.5*M*M*(GAM 


END 






THIS FUNCTION DETERMINES STATIC TEMPERATURE GIVEN MACH, GAMMA *** 
AND STAGNATION TEMPERATURE. *** 


64 
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FUNCTION FNT(TT,GAM,M) 


REAL TT,GAM,M 


FNT - TT/(1.+.5*M*M*(GAM - 1.)) 


*** THIS FUNCTION DETERMINES GAS VELOCITY GIVEN GAMMA, GAS CONSTANT *** 


*** STATIC TEMPERATURE AND MACH. 


FUNCTION FNU(GAM,R,T,M) 


REAL GAM,R,T,M,GC 


GC - 32.174 


FNU - M*SQRT(GC*R*GAM*T) 



THIS FUNCTION DETERMINES THE RATIO A/A* GIVEN MACH AND GAMMA. *** 


FUNCTION FNA(GAM.M) 


REAL GAM,M 


FNA - (1.0/M)*((2.0/(GAM+l.0))*(1.0+.5*M*M*(GAM 

+1.0)/(2.0*(GAM - 1.0))) 


1.0)))**((GAM 


*** THIS FUNCTION IS THE FUNCTION OF GAMMA AND R IN THE CONTINUITY *** 


*** EQUATION. 


FUNCTION FNGAM(GAM.R) 


REAL GAM,R,GC 


GC - 32.174 


FNGAM = SQRT((GAM*GC/R)*(2.0/(GAM+1.0))**((GAM + l.)/(GAM - 1.))) 


*** THIS FUNCTION DETERMINES THE STATIC PRESSURE GIVEN GAMMA, MACH ** 


*** AND STAGNATION PRESSURE. 


































FUNCTION FNP(M,PT,GAM) 

REAL M,GAM,PT 

FNP - PT/(1.0 + .5*M*M*(GAM - 1.0))**(GAM/(GAM 


SAMPLE OUTPUT: 
run mixer fortran 

VS FORTRAN COMPILER ENTERED. 14:12:36 

**MIXER** END OF COMPILATION 1 ****** 

**FNTT** END OF COMPILATION 2 ****** 

**FNT** END OF COMPILATION 3 ****** 

**FNU** END OF COMPILATION 4 ****** 

**FNA** END OF COMPILATION 5 ****** 

**FNGAM** END OF COMPILATION 6 ****** 

**FNP** END OF COMPILATION 7 ****** 

**FNPT** END OF COMPILATION 8 ****** 

VS FORTRAN COMPILER EXITED. 14:12:54 


EXECUTION BEGINS... 

INPUT R4,GAMMA4,TT4TH 
? 

54.415,1.258,5292 
INPUT CP3A,CP3B,CP4 
? 

.248,.332,.341 

INPUT MDOT3A.MDOT3B,TT3A,TT3B 
? 

.5,.5661,3198,4849.5 
INPUT P3A,P3B,A3A,A3B 
? 

16.4,16.4,.517.1.103 
INPUT U3A.U3B 
? 

4666.61,5594.2 
INPUT M4 
? 

1.423 


TT3BAR - 

TT4TH - 

TT4 

T4 

U4 

PT4 

P4 

V4 

U4-V4 - 

MDOT4* - 
MDOT4 - 


4193.020 

5292.000 

5182.102 

4108.812 

4280.703 

106.573 

34.372 

4280.520 

0.184 

1.0660 

1.0661 


CONTINUE? (Y/N) 
n 

R; T-0.58/1.46 14:13:56 



















APPENDIX D - SFRJ Properties 
TABLE 4 

PROPERTIES THROUGH SFRJ 


t 

.50 

.70 

1.00 

1.30 

1.80 

^ AT 

.82 

.75 

. 65 

.75 

.90 

m 0 

1.00 

1.00 

1.00 

1.00 

1.00 

“f 

.04 

. 05 

.07 

.10 

.13 

r 2 

3198.00 

3198.00 

3198.00 

3198.00 

3198.00 


4459.70 

4611.40 

4598.70 

4838.00 

4849.50 

m 4 

.27 

.28 

.29 

.32 

.34 

t 4 

4416.80 

4564.10 

4548.40 

4776.40 

4773.90 

Po 

637.87 

637.87 

637.87 

637.87 

637.87 

p 2 

146.73 

146.73 

146.73 

146.73 

146.73 

p 4 

138.79 

138.32 

137.89 

136.93 

135.70 

Po 

.40 

.40 

.40 

.40 

.40 

p 2 

142.89 

142.89 

142.89 

142.89 

142.89 

p 4 

132.56 

131.58 

130.68 

128.65 

126.07 

u 0 

5805.80 

5805.80 

5805.80 

5805.80 

5805.80 

u 2 

538.37 

538.37 

538.37 

538.37 

538.37 

u 4 

835.77 

889.56 

930.98 

1043.38 

1170.58 

u e 

7015.31 

7201.00 

7287.61 

7566.49 

7727.80 

F 

45.59 

54.86 

62.68 

77.16 

91.47 

S 

8.0492e-4 

9.3650e-4 

1.1710e-3 

1.2366e-3 

1.4444e-3 


Remarks: 

1) M 0 = 6.0; h = 80,000 ft 

2) r/ AT - fn( $ ) 

3) All quantites in standard British units 






PROPERTIES THROUGH SFRJ 


$ 

. 50 

.70 

1.00 

1.30 

1.50 

7 AT 

.82 

.75 

.65 

.75 

.81 

m 0 

1.00 

1.00 

1.00 

1.00 

1.00 

Bf 

.04 

.05 

.07 

. 10 

. 11 

t 2 

3198.00 

3198.00 

3198.00 

3198.00 

3198.00 

t 4 

4513.40 

4716.80 

4736.50 

4978.20 

5007.70 

m 4 

.27 

.29 

.30 

.32 

.33 

U 

4469.80 

4668.00 

4684.10 

4913.70 

4936.80 

Po 

637.87 

637.87 

637.87 

637.87 

637.87 

►u 

K> 

637.87 

637.87 

637.87 

637.87 

637.87 

P 4 

603.09 

600.82 

598.77 

594.47 

592.16 

Po 

.40 

.40 

.40 

.40 

.40 

P2 

621.18 

621.18 

621.18 

621.18 

621.18 

P 4 

575.75 

571.00 

566.72 

557.64 

552.74 

u 0 

5805.80 

5805.80 

5805.80 

5805.80 

5805.80 

u 2 

538.37 

538.37 

538.37 

538.37 

538.37 

u 4 

843.85 

905.64 

952.10 

1067.93 

1125.72 

u e 

7487.60 

7736.89 

7857.18 

8156.33 

8244.19 

F 

60.81 

72.38 

81.68 

97.25 

104.00 

S 

6.0349e-4 

7.0991e-4 

8.9859e-4 

9.8122e-4 

1.0587e-3 


Remarks: 

1) M 0 - 6.0; h - 80,000 ft 

2) ij at ■ fn( * ) 

3) All quantites in standard British units 

4) Isentropic inlet 















TABLE 6 





PROPERTIES THROUGH 

SFRJ 


$ 

.70 

.90 

1.10 

1.30 

1.80 

7 AT 

.75 

.68 

.68 

.75 

.90 

m 0 

1.00 

1.00 

1.00 

1.00 

1.00 

m f 

.05 

.07 

. 08 

.10 

.13 

t 2 

1638.00 

1638.00 

1638.00 

1638.00 

1638.00 

t 4 

3567. 4 0 

3655.10 

37-49.80 

3881.80 

3813.00 

m 4 

.37 

.39 

. 4 1 

. 44 

. -48 

fc 4 

350 4 .30 

3586.50 

3672.10 

3789.20 

3696.50 

Po 

61.3-4 

61.3-4 

61.3-4 

61.3-4 

61.3-4 

?2 

31.95 

31.95 

31.95 

31.95 

31.95 

P 4 

29.26 

29.11 

28.91 

28.63 

28.18 

Po 

. 4 0 4 

. 4 0 4 

. 4 0 4 

.-40-4 

.-40-4 

P2 

31.09 

31.09 

31.09 

31.09 

31.09 

P 4 

26.85 

26.51 

26.07 

25. 4 5 

2-4. -42 

u O 

3870.55 

3870.55 

3870.55 

3870.55 

3870.55 

u 2 

390.90 

390.90 

390.90 

390.90 

390.90 

u 4 

1018.35 

1078.20 

1156.89 

1271.23 

1 4 37.11 

u e 

5785.5 4 

590 4 .23 

6030.25 

6188.78 

6272.16 

F 

68.76 

75.33 

82.26 

90. 4 1 

100.-40 

S 

7. 4 725e- 4 

8.7692e- 4 

9.815 4 e- 4 1, 

. 055 4 e-3 1. 

3159e-3 

Remarks: 






1) M 0 - 4 . 

, 0; h - 80, 

000 ft 




2) n&j = f n ( 4> ) 





3) All quantites in 

standard British units 


















APPENDIX E - DMRJ Properties 
TABLE 7 

PROPERTIES THROUGH DMRJ (50/50) 


*SFRJ 

1.20 

1.50 

1.80 

2.20 

2.30 

$tot 

.60 

.75 

.90 

1.10 

1.15 

m 3A 

. 50 

.50 

.50 

.50 

.50 

m 3B 

. 5 4 

.56 

.57 

.58 

.58 

t 3 a 

3198.00 

3198.00 

3198.00 

3198.00 

3198.00 

T 3B 

4 77 4 .50 

4 903. 4 0 

4 8 4 9.50 

<4612.80 

4 53 4 .3 0 

t 4 

4 856. 4 2 

5056.25 

5182.10 

5258.11 

5263.50 

M 3B 

2.06 

2.05 

2.05 

2 . 0 4 

2.0-4 

m 4 

1.5 4 

1. 4 9 

1. 4 2 

1.30 

1.27 

fc 3B 

2976.70 

303 4 .30 

2922.50 

2766.30 

2706.50 

t 4 

3701.93 

3927.17 

4 108.81 

4 312.30 

4 359.3 4 

p o 

637.87 

637.87 

637.87 

637.87 

637.87 

P 2A 

281.96 

281.96 

281.96 

281.96 

281.96 

P 2B 

1 4 6.73 

1 4 6.73 

1 4 6.73 

1 4 6.73 

1 4 6.73 

P 3A 

281.96 

281.96 

281.96 

281.96 

281.96 

P 3B 

138.01 

136.36 

135.72 

135.17 

135.09 

p 4 

113.06 

109.98 

106.57 

102.67 

101.72 

Po 

. 4 0 

. 4 0 

. 4 0 

. 4 0 

.<40 

P3A 

16. 4 0 

16. 4 0 

16. 4 0 

16. 4 0 

16.<40 

P3B 

16. 4 0 

16. 4 0 

16. 4 0 

16.<40 

16. 4 0 

P 4 

30.70 

32.20 

3 4 .37 

39.0-4 

4 0.58 

U 3A 

4 666.61 

4 666.61 

4 666.61 

4 666.61 

<4 666.61 

U 3B 

5391.56 

5539.50 

559 4 .2 0 

5568.70 

55<48.60 

u 4 

4 365.81 

4 630.70 

4 280.70 

<4052.73 

3971.72 

u e 

7295. 4 8 

7 4 88.00 

7620.32 

7730.2 4 

77-48.30 

u 0 

5805.80 

5805.80 

5805.80 

5805.80 

5805.80 

F 

56.28 

65.11 

72.05 

79.20 

80.70 

S 

7.8183e- 4 

8. 4 626e- 4 

9.1738e- 4 1 

.019<4e-3 1 

. 0<458e-3 


Remarks: 

1) M 0 = 6.0; h = 80,000 ft 







PROPERTIES THROUGH DMRJ (50/50) 


$ SFRJ 

1.20 

1.50 

1.80 

2.20 

2.30 

^TOT 

.60 

.75 

.90 

1.10 

1.15 

m 3A 

.50 

. 50 

.50 

.50 

.50 

m 3B 

.54 

.56 

.57 

. 58 

.58 

t 3 a 

3198.00 

3198.00 

3198.00 

3198.00 

3198.00 

t 3B 

4774.50 

4903.40 

4849.50 

4612.80 

4534.30 

t 4 

4611.70 

4774.99 

4852.41 

4859.43 

4848.02 

m 3B 

2.06 

2.05 

2.05 

2.04 

2.04 

m 4 

1.69 

1.66 

1.62 

1.55 

1.53 


2976.70 

3034.30 

2922.50 

2766.30 

2706.50 

t 4 

3352.56 

3523.64 

3628.52 

3708.57 

3721.33 

Po 

637.87 

637.87 

637.87 

637.87 

637.87 

P2A 

281.96 

281.96 

281.96 

281.96 

281.96 

P 2B 

146.73 

146.73 

146.73 

146.73 

146.73 

P 3 A 

281.96 

281.96 

281.96 

281.96 

281.96 

p 3 B 

138.01 

136.36 

135.72 

135.17 

135.09 

P 4 

123.12 

120.22 

117.27 

113.69 

112.69 

Po 

.40 

.40 

.40 

.40 

.40 

P 3 A 

16.40 

16.40 

16.40 

16.40 

16.40 

P 3 B 

16.40 

16.40 

16.40 

16.40 

16.40 

P 4 

26.63 

27.44 

28.43 

30.44 

31.03 

U 3A 

4666.61 

4666.61 

4666.61 

4666.61 

4666.61 

U 3 B 

5391.56 

5539.50 

5594.20 

5568.70 

5548.60 

u 4 

4559.37 

4590.77 

4571.16 

4470.29 

4433.61 

u e 

7135.60 

7305.35 

7405.50 

7465.70 

7470.70 

u 0 

5805.80 

5805.80 

5805.80 

5805.80 

5805.80 

F 

51.09 

59.12 

64.93 

70.33 

71.34 

S 

8.6123e-4 9, 

. 3203e-4 

1.0179e-3 1. 

. 1481e-3 1. 

. 1830e-3 


Remarks 



















PROPERTIES THROUGH DMRJ (50/50) 


4 SFRJ 

1.20 

1.50 

1.80 

2.20 

2.30 

4 TOT 

.60 

.75 

.90 

1.10 

1.15 

m 3A 

.50 

.50 

.50 

.50 

.50 

m 3B 

• 5 4 

.56 

.57 

.58 

.58 

t 3 a 

3198.00 

3198.00 

3198.00 

3198.00 

3198.00 

T 3B 

4 77 4 .50 

4 903. 4 0 

4 8 4 9.50 

4 612.80 

4 53 4 .30 

t 4 

4 530.13 

4 681.2 4 

4 7 4 2.51 

4 726.53 

4 709.53 

m 3B 

2.06 

2.05 

2.05 

2.04 

2.04 

m 4 

1 . 7 4 

1.71 

1.68 

1.63 

1.62 

^B 

2976.70 

303 4 .30 

2922.50 

2766.30 

2706.50 

t 4 

3239.07 

3392.20 

3 4 73.37 

3518.96 

3522.78 

p 0 

637.87 

637.87 

637.87 

637.87 

637.87 

P 2A 

281.96 

281.96 

281.96 

281.96 

281.96 

P 2B 

1 4 6.73 

1 4 6.73 

1 4 6.73 

1 4 6.73 

146.73 

P 3A 

281.96 

281.96 

281.96 

281.96 

281.96 

P 3B 

138.01 

136.36 

135.72 

135.17 

135.09 

p 4 

127.22 

12 4 .57 

122.00 

118.82 

117.90 

PO 

. 4 0 

. 4 0 

. 4 0 

. 4 0 

.40 

P3A 

16. 4 0 

16. 4 0 

16. 4 0 

16. 4 0 

16.40 

P3B 

16. 4 0 

16. 4 0 

16. 4 0 

16. 4 0 

16.40 

P 4 

25.-40 

26.03 

26.72 

28.19 

28.62 

U 3A 

4 666.61 

4 666.61 

4 666.61 

4 666.61 

4666.61 

U 3B 

5391.56 

5539.50 

559 4 .20 

5568.70 

5548.60 

u 4 

4 616.78 

4 659.38 

4 65 4 .91 

4 579.11 

4550.24 

u e 

7082.10 

72 44 . 4 0 

7333.80 

7377.30 

7378.00 

u 0 

5805.80 

5805.80 

5805.80 

5805.80 

5805.80 

F 

4 9.35 

57.12 

62.56 

67.36 

68.22 

S 

8.9152e- 4 

9.6 4 6 4 e- 4 

1.0566e-3 

1.1987e-3 

1.2372e-3 


Remarks 











TABLE 10 

PROPERTIES THROUGH DMRJ (50/50) 


*SFRJ 

1.20 

1.50 

1.80 

2.20 

2.30 

*TOT 

.60 

.75 

.90 

1.10 

1.15 

m 3A 

.50 

.50 

.50 

.50 

.50 

m 3B 

.54 

.56 

.57 

.58 

.58 

t 3 a 

3198.00 

3198.00 

3198.00 

3198.00 

3198.00 

T 3B 

4774.50 

4903.40 

4849.50 

4612.80 

4534.30 

t 4 

4326.19 

4446.86 

4467.76 

4394.30 

4363.30 

M 3B 

2.06 

2.05 

2.05 

2.04 

2.04 

m 4 

1.88 

1.86 

1.86 

1.84 

1.84 

t^B 

2976.70 

3034.30 

2922.50 

2766.30 

2706.50 

fc 4 

2957.48 

3068.96 

3094.25 

3058.51 

3040.95 

p 0 

637.87 

637.87 

637.87 

637.87 

637.87 

P 2A 

281.96 

281.96 

281.96 

281.96 

281.96 

P 2B 

146.73 

146.73 

146.73 

146.73 

146.73 

P 3A 

281.96 

281.96 

281.96 

281.96 

281.96 

P 3B 

138.01 

136.36 

135.72 

135.17 

135.09 

P 4 

139.95 

138.17 

137.20 

136.36 

136.12 

Po 

.40 

.40 

.40 

.40 

.40 

P3A 

16.40 

16.40 

16.40 

16.40 

16.40 

P3B 

16.40 

16.40 

16.40 

16.40 

16.40 

P4 

22.53 

22.78 

22.88 

23.30 

23.41 

U 3A 

4666.61 

4666.61 

4666.61 

4666.61 

4666.61 

U 3B 

5391.56 

5539.50 

5594.20 

5568.70 

5548.60 

u 4 

4753.61 

4817.32 

4842.54 

4816.07 

4803.16 

u e 

6948.42 

7091.72 

7153.93 

7155.51 

7145.50 

u 0 

5805.80 

5805.80 

5805.80 

5805.80 

5805.80 

F 

45.02 

52.11 

56.60 

59.91 

60.38 

S 

9.7743e-4 1. 

. 0573e-3 

1.1679e-3 1, 

. 3478e-3 1, 

, 3977e-3 


Remarks: 


1) M 0 - 6.0; h = 80,000 ft 

2) “ 0.25 







$ SFRJ 

TABLE 11 

PROPERTIES THROUGH DMRJ 

1.10 1.30 1.50 

(80/20) 

1.80 

2.20 

4 TOT 

.88 

1.0 4 

1.20 

1.44 

1.76 

m 3A 

.20 

.20 

.20 

.20 

.20 

m 3B 

.86 

.88 

.89 

.91 

.93 

t 3 a 

3198.00 

3198.00 

3198.00 

3198.00 

3198.00 

T 3B 

4 693. 4 0 

4 838.00 

4 903. 4 0 

4849.50 

4612.80 

t 4 

519 4 .71 

5293.01 

5330.93 

5271.16 

5007.53 

M 3B 

2.057 

2.05 

2.05 

2.05 

2.04 

m 4 

1.56 

1.58 

1.60 

1.61 

1.62 

t 3B 

2928. 4 0 

3010.20 

303 4 .30 

2963.10 

2753.00 

t 4 

3958. 4 0 

4 001.22 

4009.98 

3936.22 

3714.24 

p o 

637.87 

637.87 

637.87 

637.87 

637.87 

P 2A 

281.96 

281.96 

281.96 

281.96 

281.96 

P 2B 

1 4 6.73 

1 4 6.73 

146.73 

146.73 

146.73 

P 3A 

281.96 

281.96 

281.96 

281.96 

281.96 

P 3B 

137.56 

136.92 

136.35 

135.71 

135.16 

p 4 

103.60 

10 4 .79 

105.61 

106.43 

107.50 

PO , 

. 4 0 

.40 

.40 

.40 

.40 

P3A 

16. 4 0 

16. 4 0 

16.40 

16.40 

16.40 

P3B 

16. 4 0 

16. 4 0 

16.40 

16.40 

16.40 

P 4 

27.53 

26.78 

26.35 

26.07 

26.04 

U 3A 

4 666.61 

4 666.61 

4666.61 

4666.61 

4666.61 

U 3B 

5323. 4 7 

5 4 50.17 

5539.40 

5594.10 

5562.50 

U 4 

4 590.69 

4723. 4 3 

4812.61 

4873.04 

4859.37 

u e 

7607.70 

7738.90 

7823.28 

7856.59 

7736.12 

u 0 

5805.80 

5805.80 

5805.80 

5805.80 

5805.80 

F 

71.28 

78. 44 

84.13 

89.55 

91.06 

S 9 

. 0629e-4 9 

.7278e- 4 

1.0472e-3 

1.1803e-3 

1.4188e-3 

Remarks: 

1) M 0 

- 6.0; h 

- 80,000 

ft 


















TABLE 12 




PROPERTIES THROUGH DMRJ 

(80/20) 


*SFRJ 

1.10 

1.30 

1.50 

1.80 

2.20 

$tot 

.88 

1.0 4 

1.20 

1. 44 

1.76 

m 3A 

.20 

.20 

.20 

.20 

.20 

m 3B 

.86 

.88 

.89 

.91 

.93 

t 3 a 

3198.00 

3198.00 

3198.00 

3198.00 

3198.00 

T 3B 

4 693. 4 0 

4 838.00 

4 903. 4 0 

4 8 4 9.50 

4 612.80 

t 4 

4 952.03 

5060.26 

5105.70 

5052.66 

4 811.1 4 

M 3B 

2.06 

2.05 

2.05 

2.05 

2.0 4 

m 4 

1.69 

1.71 

1.72 

1.73 

1.73 

fc 3B 

2928. 4 0 

3010.20 

303 4 .30 

2963.10 

2753.00 

t 4 

361 4 .1 4 

3673.11 

3693.0 4 

3628. 4 2 

3 4 35.36 

p o 

637.87 

637.87 

637.87 

637.87 

637.87 

*rJ 

to 

> 

281.96 

281.96 

281.96 

281.96 

281.96 

P 2B 

1 4 6.73 

1 4 6.73 

1 4 6.73 

1 4 6.73 

1 4 6.73 

P 3A 

281.96 

281.96 

281.96 

281.96 

281.96 

P 3B 

137.56 

136.92 

136.35 

135.71 

135.16 

p 4 

112.21 

113.15 

113.8 4 

11 4 .66 

115. 4 8 

PO 

. 4 0 

. 4 0 

. 4 0 

. 4 0 

. 4 0 

P3A 

16. 4 0 

16. 4 0 

16. 4 0 

16. 4 0 

16. 4 0 

p 3B 

16. 4 0 

16. 4 0 

16. 4 0 

16. 4 0 

16. 4 0 

P 4 

2 4 .16 

23.73 

23. 4 6 

23.27 

23.36 

U 3A 

4 666.61 

4 666.61 

4 666.61 

4 666.61 

4 666.61 

U 3B 

5323. 4 7 

5 4 50.17 

5539. 4 0 

559 4 .10 

5562.50 

u 4 

4 775.56 

4 89 4 .66 

4 976.88 

5033.38 

-5011.57 

u e 

7 4 5 4 .65 

7593.05 

7682.08 

7717.50 

7606.55 

u o 

5805.80 

5805.80 

5805.80 

5805.80 

5805.80 

F 

66.22 

73.56 

79.35 

8 4 .77 

86.51 

S 9, 

. 7560e- 4 1. 

0373e-3 

1.1102e-3 

1.2 4 69e-3 1, 

. 4 93 4 e-3 

Remarks: 






1) H 0 

= 6.0; h = 

» 80,000 

ft 



2 ) 7 at 

= 0.6 


75 










VV,< 




a 


’S’*! 


TABLE 13 

PROPERTIES THROUGH DMRJ (80/20) 


*SFRJ 

1.10 

1.30 

1.50 

1.80 

2.20 

$tot 

.88 

1.0 4 

1.20 

1.44 

1.76 

m 3A 

.20 

.20 

.20 

.20 

.20 

m 3B 

.86 

.88 

.89 

.91 

.93 

t 3 a 

3198.00 

3198.00 

3198.00 

3198.00 

3198.00 

t 3 b 

4 693. 4 0 

■4838.00 

■4903. 4 0 

4849.50 

4612.80 

t 4 

4 871.13 

4 982.67 

5030.63 

4979.82 

4745.50 

M 3B 

2.06 

2.05 

2.05 

2.05 

2.04 

m 4 

1 . 7 4 

1.76 

1.77 

1.77 

1.77 

t 3B 

2928. *40 

3010.20 

303 4 .30 

2963.10 

2753.00 

t 4 

3501.87 

3565.87 

3588.53 

3527.28 

3344.93 

p o 

637.87 

637.87 

637.87 

637.87 

637.87 

P 2A 

281.96 

281.96 

281.96 

281.96 

281.96 

P 2B 

1 4 6.73 

1 4 6.73 

1 4 6.73 

146.73 

146.73 

P 3A 

281.96 

281.96 

281.96 

281.96 

281.96 

P 3B 

137.56 

136.92 

136.35 

135.71 

135.16 

p 4 

115.68 

116. 4 7 

117.16 

117.93 

118.50 

PO 

. 4 0 

.40 

.40 

.40 

.40 

P3A 

16. 4 0 

16. 4 0 

16.40 

16.40 

16.40 

P3B 

16. 4 0 

16. 4 0 

16.40 

16.40 

16.40 

P 4 

23.1 4 

22.79 

22.57 

22.40 

22.54 

U 3A 

■4 666.61 

4 666.61 

4666.61 

4666.61 

4666.61 

U 3B 

5323. 4 7 

5 4 50.17 

5539.40 

5594.10 

5562.50 

u 4 

4 831.23 

4 9 4 6.70 

5028.47 

5083.15 

5056.53 

u e 

7 4 03.50 

75 44 .30 

7635.10 

7671.05 

7562.83 

u 0 

5805.80 

5805.80 

5805.80 

5805.80 

5805.80 

F 

6 4 .52 

71.93 

77.76 

83.18 

84.98 

S 

1.0012e-3 1, 

. 0608e-3 

1.1329e-3 1, 

, 2708e-3 

1.5204e-3 


Remarks: 

1) M 0 ■ 6.0; h 

2) r? AT « 0.5 


80,000 ft 















PROPERTIES THROUGH DMRJ (80/20) 


$ SFRJ 

1.10 

1.30 

1.50 

1.80 

2.20 

4 TOT 

.88 

1.04 

1.20 

1.44 

1.76 

m 3A 

.20 

.20 

.20 

.20 

.20 

m 3B 

.86 

.88 

.89 

.91 

.93 

t 3 a 

3198.00 

3198.00 

3198.00 

3198.00 

3198.00 

T 3B 

4693.40 

4838.00 

4903.40 

4849.50 

4612.80 

t 4 

4668.90 

4788.71 

4842.95 

4797.74 

4581.75 

m 3B 

2.06 

2.05 

2.05 

2.05 

2.04 

m 4 

1.86 

1.87 

1.88 

1.88 

1.88 

fc 3B 

2928.40 

3010.20 

3034.30 

2963.10 

2753.00 

U 

3223.91 

3300.05 

3331.88 

3277.16 

3118.25 

P 0 

637.87 

637.87 

637.87 

637.87 

637.87 

P 2A 

281.96 

281.96 

281.96 

281.96 

281.96 

P 2B 

146.73 

146.73 

146.73 

146.73 

146.73 

P 3A 

281.96 

281.96 

281.96 

281.96 

281.96 

P 3B 

137.56 

136.92 

136.35 

135.71 

135.16 

P 4 

126.18 

126.41 

126.77 

127.55 

127.62 

Po 

.40 

.40 

.40 

.40 

.40 

P3A 

16.40 

16.40 

16.40 

16.40 

16.40 

P3B 

16.40 

16.40 

16.40 

16.40 

16.40 

P4 

20.74 

20.58 

20.47 

20.34 

20.55 

U 3A 

4666.61 

4666.61 

4666.61 

4666.61 

4666.61 

U 3B 

5323.47 

5450.17 

5539.40 

5594.10 

5562.50 

u 4 

4963.02 

5070.58 

5147.29 

5200.83 

5168.88 

u e 

7275.70 

7422.45 

7517.00 

7554.90 

7454.49 

u 0 

5805.80 

5805.80 

5805.80 

5805.80 

5805.80 

F 

60.29 

67.85 

73.77 

79.18 

81.18 

S 

1.0714e-3 1, 

. 1246e-3 

1.1943e-3 1, 

,3349e-3 1. 

. 5916e-3 


Remarks 






J 


Si 




TABLE 15 




PROPERTIES THROUGH DMRJ 

(80/20) 


*SFRJ 

1.20 

1.50 

1.80 

2.00 

2.20 

$tot 

.96 

1.20 

1.44 

1.60 

1.76 

m 3A 

.20 

.20 

.20 

.20 

.20 

m 3B 

.87 

.89 

.91 

.92 

.93 

t 3 a 

1638.00 

1638.00 

1638.00 

1638.00 

1638.00 

T 3B 

3837.10 

3882.50 

3813.00 

3728.00 

3620.60 

t 4 

4540.29 

4602.13 

4394.83 

4215.56 

4027.05 

M 3B 

2.26 

2.25 

2.25 

2.25 

2.25 

m 4 

1.32 

1.33 

1.41 

1.47 

1.53 

t 3B 

2236.70 

2222.40 

2140.00 

2065.10 

1979.60 

t 4 

3718.25 

3754.16 

3481.94 

3262.25 

3051.32 

p o 

61.34 

61.34 

61.34 

61.34 

61.34 

P 2A 

41.77 

41.77 

41.77 

41.77 

41.77 

P 2B 

31.95 

31.95 

31.95 

31.95 

31.95 

P 3A 

41.77 

41.77 

41.77 

41.77 

41.77 

P 3B 

28.76 

28.42 

28.18 

28.05 

27.95 

P 4 

15.94 

16.05 

16.61 

17.09 

17.56 

1 Po 

.40 

.40 

.40 

.40 

.40 

P3A 

2.44 

2.44 

2.44 

2.44 

2.44 

P3B 

2.44 

2.44 

2.44 

2.44 

2.44 

p 4 

5.93 

5.89 

5.43 

5.10 

4.85 

! U 3A 

3316.50 

3316.50 

3316.50 

3316.50 

3316.50 

U 3B 

5080.44 

5193.07 

5222.49 

5216.18 

5189.83 

U 4 

3744.51 

3837.09 

4002.05 

4100.34 

4160.59 

u e 

6297.70 

6405.80 

6359.22 

6284.75 

6205.44 

; u o 

3870.55 

3870.55 

3870.55 

3870.55 

3870.55 

F 

89.24 

96.34 

98.24 

97.97 

97.49 

! 

. 9002e-4 9 

.1448e-4 

1.0759e-3 

1.1983e-3 

1.3253e-3 

t 

] Remarks: 






! 1) M 0 

* 4.0; h > 

= 80,000 

ft 



; 2 ) 7 at 

! 

! 

> 

\ 

* 

l 

t - . ... • 

= 0.9 


78 

, rm jv , 

k_*v _*v jv juci i. 

V -V JV .V 










TABLE 16 

PROPERTIES THROUGH DMRJ (80/20) 


$ SFRJ 

1.20 

1.50 

1.80 

2.00 

2.20 

$tot 

.96 

1.20 

1.44 

1.60 

1.76 

m 3A 

.20 

.20 

.20 

.20 

.20 

m 3B 

.87 

.89 

.91 

.92 

.93 

t 3 a 

1638.00 

1638.00 

1638.00 

1638.00 

1638.00 

T 3B 

3837.10 

3882.50 

3813.00 

3728.00 

3620.60 

t 4 

4197.95 

4254.53 

4098.51 

3956.26 

3801.00 

M 3B 

2.26 

2.25 

2.25 

2.25 

2.25 

m 4 

1.56 

1.57 

1.62 

1.66 

1.70 

^B 

2236.70 

2222.40 

2140.00 

2065.10 

1979.60 

U 

3206.91 

3236.18 

3046.01 

2884.58 

2723.74 

p 0 

61.34 

61.34 

61.34 

61.34 

61.34 

P 2A 

41.77 

41.77 

41.77 

41.77 

41.77 

P 2B 

31.95 

31.95 

31.95 

31.95 

31.95 

P 3A 

41.77 

41.77 

41.77 

41.77 

41.77 

P 3B 

28.76 

28.42 

28.18 

28.05 

27.95 

P 4 

17.67 

17.81 

18.40 

18.87 

19.33 

Po 

.40 

.40 

.40 

.40 

.40 

P3A 

2.44 

2.44 

2.44 

2.44 

2.44 

P3B 

2.44 

2.44 

2.44 

2.44 

2.44 

P4 

4.66 

4.63 

4.42 

4.25 

4.12 

U 3A 

3316.50 

3316.50 

3316.50 

3316.50 

3316.50 

U 3B 

5080.44 

5193.07 

5222.49 

5216.18 

5189.83 

U 4 

4111.46 

4204.95 

4297.20 

4347.46 

4371.70 

u e 

6110.30 

6214.80 

6194.14 

6138.50 

6075.57 

u 0 

3870.55 

3870.55 

3870.55 

3870.55 

3870.55 

F 

83.00 

89.88 

92.57 

92.89 

92.93 

S 

8.4937e-4 9 

.8020e-4 

1.1419e-3 1, 

, 2639e-3 1. 

. 3903e-3 

Remarks: 






1) 

M 0 = 4.0; h 

= 80,000 

ft 




2) *?at “ 0.6 























TABLE 17 




PROPERTIES THROUGH DMRJ 

(75/25) 


S’SFRJ 

1.10 

1.30 

1.50 

1.80 

2.20 

$tot 

.825 

.975 

1.125 

1.35 

1.65 

m 3A 

.25 

.25 

.25 

.25 

.25 

m 3B 

.81 

. 82 

.83 

.85 

.87 

t 3 a 

3198.00 

3198.00 

3198.00 

3198.00 

3198.00 

T 3B 

4693.40 

4838.00 

4903.40 

4849.50 

4612.80 

t 4 

5144.86 

5257.99 

5315.18 

5300.77 

5101.19 

M 3B 

2.06 

2.05 

2.05 

2.05 

2.04 

m 4 

1.55 

1.56 

1.57 

1.55 

1.55 

fc 3B 

2928.80 

3010.60 

3034.90 

2963.70 

2767.00 

t 4 

3930.61 

4000.49 

4034.09 

4035.50 

3864.15 

p o 

637.87 

637.87 

637.87 

637.87 

637.87 

P 2A 

281.96 

281.96 

281.96 

281.96 

281.96 

P 2B 

146.74 

146.74 

146.74 

146.74 

146.74 

P 3A 

281.96 

281.96 

281.96 

281.96 

281.96 

P 3B 

137.48 

136.83 

136.25 

135.58 

135.02 

P 4 

104.58 

105.17 

105.50 

104.63 

104.93 

PO 

.40 

.40 

.40 

.40 

.40 

P3A 

16.40 

16.40 

16.40 

16.40 

16.40 

P3B 

16.40 

16.40 

16.40 

16.40 

16.40 

P4 

28.20 

27.74 

27.50 

27.91 

28.04 

U 3A 

4665.61 

4666.61 

4666.61 

4666.61 

4666.61 

U 3B 

5322.90 

5449.47 

5538.60 

5593.17 

5567.56 

u 4 

4536.71 

4647.34 

4722.08 

4733.79 

4721.40 


7552.75 

7696.37 

7793.20 

7849.98 

7826.360 

U 0 

5805.80 

5805.80 

5805.80 

5805.80 

5805.80 

F 

68.52 

75.89 

81.78 

87.71 

92.26 

S 8, 

. 8438e-4 9 

.4349e-4 

1.0101e-3 

1.1298e-3 

1.3126e-3 

Remarks: 






1) M 0 

= 6.0; h 

= 80,000 

ft 















TABLE 18 

PROPERTIES THROUGH DMRJ (60/40) 


$ SFRJ 


1.30 

1.50 

1.80 

2.20 

$ T ot 


.78 

.90 

1.08 

1.32 

m 3A 


.40 

.40 

.40 

.40 

m 3B 


. 66 

.67 

. 68 

.70 

t 3 a 


3198.00 

3198.00 

3198.00 

3198.00 

T 3B 


4838.00 

4903.40 

4849.50 

4612.80 

t 4 


5095.59 

5199.88 

5281.00 

5271.89 

m 3B 


2.05 

2.05 

2.05 

2.04 

m 4 


1.52 

1.51 

1.46 

1.37 

^B 


3010.20 

3034.30 

2963.10 

2766.30 

t 4 


3917.36 

4024.09 

4414.48 

4239.85 

P 0 


637.87 

637.87 

637.87 

637.87 

P 2A 


281.96 

281.96 

281.96 

281.96 

P 2B 


146.74 

146.74 

146.74 

146.74 

P 3A 


281.96 

281.96 

281.96 

281.96 

P 3B 


136.92 

136.35 

135.71 

135.16 

P 4 


108.15 

106.97 

104.78 

101.54 

Po 


.40 

.40 

.40 

.40 

P3A 


16.40 

16.40 

16.40 

16.40 

P3B 


16.40 

16.40 

16.40 

16.40 

P4 


30.12 

30.65 

32.14 

35.21 

U 3A 


4666.61 

4666.61 

4666.61 

4666.61 

U 3B 


5450.16 

5539.41 

5594.11 

5568.60 

u 4 


4458.84 

4480.46 

4438.38 

4274.51 

u e 


7518.44 

7634.66 

7742.18 

7800.58 

u 0 


5805.80 

5805.80 

5805.80 

5805.80 

F 


66.62 

72.53 

79.27 

85.52 

S 

8, 

. 6010e-4 9 

.1137e-4 

1.0004e-3 

1.1343e-3 

$ 

Remarks: 






1) 

Mq 

= 6.0; h 

= 80,000 

ft 



2.30 

1.38 

.40 

.70 

3198.00 

4534.30 

5249.18 

2.04 

1.35 

2706.50 

4240.61 

637.87 
281.96 
146.74 
281.96 
135.08 

100.88 

.40 

16.40 

16.40 

35.98 

4666.61 
5548.52 
4227.20 
7815.95 
5805.80 

87.09 
1.1632e-3 


2) V AT ~ 0*9 




mm 














PROPERTIES THROUGH DMRJ (50/50) 


*SFRJ 

1.20 

1.50 

1.80 

2.00 

2.20 

$tot 

.60 

.75 

.90 

1.00 

1.10 

m 3A 

.50 

.50 

.50 

.50 

.50 

m 3B 

.54 

.56 

.57 

.57 

.58 

T 3A 

3198.00 

3198.00 

3198.00 

3198.00 

3198.00 

T 3B 

4923.00 

5007.70 

4906.30 

4785.70 

4638.40 

t 4 

4959.27 

5206.12 

5364,06 

5425.81 

5454.29 

M 3B 

2.44 

2.44 

2.44 

2.44 

2.44 

m 4 

1.74 

1.66 

1.57 

1.51 

1.45 

fc 3B 

2737.00 

2749.70 

2642.50 

2544.80 

2435.30 

fc 4 

3558.25 

3849.85 

4087.70 

4202.44 

4307.33 

p o 

637.87 

637.87 

637.87 

637.87 

637.87 

P 2A 

637.87 

637.87 

637.87 

637.87 

637.87 

P 2B 

637.87 

637.87 

637.87 

637.87 

637.87 

P 3A 

637.87 

637.87 

637.87 

637.87 

637.87 

P 3B 

596.23 

592.67 

590.18 

588.97 

588.06 

p 4 

343.50 

327.83 

312.65 

303.03 

297.10 

Po 

.40 

.40 

.40 

.40 

.40 

P3A 

37.08 

37.08 

37.08 

37.08 

37.08 

P3B 

37.08 

37.08 

37.08 

37.08 

37.08 

P4 

69.08 

74.57 

81.74 

85.83 

92.62 

U 3A 

4666.61 

4666.61 

4666.61 

4666.61 

4666.61 

U 3B 

6072.34 

6215.17 

6248.42 

6234.66 

6198.93 

u 4 

4813.48 

4784.97 

4677.20 

4596.10 

4470.48 

u e 

7760.03 

8000.26 

8153.93 

8236.10 

8276.85 

u 0 

5805.80 

5805.80 

5805.80 

5805.80 

5805.80 

F 

71.38 

81.91 

89.73 

94.33 

97.56 

S 

6.178 6e-4 

6.7271e-4 

7.3662e-4 7, 

ft 

,7816e-4 8. 

. 2716e-4 

Remarks 






1) 

M 0 = 6.0; h 

= 80,000 

ft 



2) 

7 AT = 0 • 9 





3) 

Isentropic 

inlet 
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